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Abstract. A natural approach to measure the time of arrival of an atom at a spatial 
region is to illuminate this region with a laser and detect the first fluorescence photons 
produced by the excitation of the atom and subsequent decay. We investigate the 
actual physical content of such a measurement in terms of atomic dynamical variables, 
taking into account the finite width of the laser beam. Different operation regimes 
are identified, in particular the ones in which the quantum current density may be 
obtained. 
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1. Introduction 

In spite of the emphasis on "observables" and "measurements" in the standard 
formulation of quantum mechanics, one of the fundamental difficulties of the theory 
is to establish the relation between "observables" represented by operators that act on 
the variables of the microscopic system in isolation, and actual measurements. There is 
no universal prescription to construct an apparatus for a given operator and, in the other 
direction, identifying from the macroscopic experiment the corresponding microscopic 
operator may also be difficult. These problems have been particularly accute, or have 
been perceived to be so, for time observables. According to an argument due to Pauli 
PP, there is no self-adjoint time operator conjugate to a semibounded Hamiltonian 
(see however 0), and this limitation has hindered the research on the theoretical 
understanding and description of time quantities measured in the laboratories. The 
trend is slowly changing though During the last two decades, many works have 
been devoted to find consistent theories, in terms of operators or otherwise, about 
characteristic times for tunneling, quantum jumps, decay, or arrival jUIH|. This effort 
has been mostly focused on establishing fundamental theories for the isolated system, 
typically a structureless particle. In these theories the measuring apparatus is totally 
absent, at least in an explicit form, or at best modelled very schematically jSHH]- We have 
advocated the need for an analysis of more realistic models [H Ej , not only to determine 
how and if the ideally defined quantities can be measured but also to ascertain what 
usual operational methods are "really" measuring in the language of system variables. 
In a previous article [7j we have modelled the measurement by fluorescence of the time 
of arrival (TOA) of an atom at a given location. An important finding was that, in 
the limit of weak driving, the quantum current density J may be obtained from the 
distribution of first photon detections II. This is satisfactory because, in the classical 
limit, J is the correct TOA distribution for the initial states considered (with momenta 
fully directed towards the laser); but it is also surprising because J may be negative in 
quantum mechanics, even for states formed entirely by positive momentum components 
jHJ. By contrast, the positive and axiomatic TOA distribution of Kijowski [5] was not 
obtained in any limit of the parameters involved, except when it is sufficiently close to 
J. A simplifying assumption of the model was that the laser-illuminated region had a 
semi-infinite width, from x = to oo. This is a good approximation at very low kinetic 
energies, because the atom is either detected before abandoning the finite illuminated 
region, or reflected by the laser; but it fails otherwise. In this paper we remove that 
assumption and consider the more realistic case of a finite laser beam width, L. 

There were two basic difficulties with a direct "good TOA measurement" in the 
semi-infinite model: the atomic reflection with no photon emission, due to a strong laser 
field, and the detection delay for a weak one. Although we could find and characterize 
operation regimes where the two effects were negligible in practice, the best fit with 
the free-motion current density was obtained for weak driving and for a short lifetime 
I/7 of the excited state (7 is Einstein's coefficient); this avoids reflection but induces 
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delay because, in spite of the short life time, the laser pumping is very inefficient. 
Nevertheless, the flux J for the freely moving atoms may be obtained in principle from 
the experimental signal n, since, in that limit, IT = J * W, where W is the (known) 
distribution of photon detection times corresponding to a laser-illuminated atom at rest. 

The applicability of the above "deconvolution" method is limited because the long 
detection delays implied require also a large laser-beam width, and the required width 
increases with the atomic velocity. In practice the laser beam is of course always finite, so 
there is a third problem aside from reflection and delay: atomic transmission through the 
laser without photon detection, which mainly affects the fast components and may also 
cause distortion in the measured photon distribution. We shall characterize this effect, 
and describe several operation regimes. Finally, we shall see how, taking advantage of 
the Rabi oscillation induced by the laser, the flux may also be obtained by deconvolution 
with a finite width laser at least in two cases: for strong driving and 7 — *> 0, i.e., in a 
limit which is quite the opposite from the one used for the semi-infinite laser, and 
when the momentum width of the atomic wave packet is small compared to its average 
momentum. 

2. The model 

The setting of the modelled experiment and the fundamental theory are described in 
Ref. so we shall only outline the basic ideas and equations here, emphasizing the 
novelties due to the finite laser-beam width. A two-level atom wave packet impinges on 
a perpendicular laser beam at resonance with the atomic transition. In the so called 
quantum jump approach ^Hj the continuous measurement of the fluorescence photons 
is simulated by a periodic projection onto no-photon or one-photon subspaces every 
St, a time interval large enough to avoid the Zeno effect, but smaller than any other 
characteristic time. The amplitude for the undetected atoms in the interaction picture 
for the internal Hamiltonian obeys, in a time scale coarser than St, and using the 
rotating wave and dipole approximations, an effective Schrodinger equation governed 
by the complex "conditional" Hamiltonian (the hat is used to distinguish momentum 
and position operators from the corresponding c-numbers) 



and f2 is the Rabi frequency. The probability, N t , of no photon detection from t Q , the 
instant when the packet is prepared far from the laser and with positive momenta, up 
to time t, is given by |T0] 





< x < L 
otherwise 



(2) 




(3) 
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and the probability density, II(i), for the first photon detection by 

n W = -§ = 7 P 2 , (4) 

where P2 is the population of the excited state. II plays the role of an "operational" 
time-of-arrival distribution. To obtain the time development under H c of a wave packet 
incident from the left we first solve the stationary equation 

H C * = E<S>, where *(x) = (^j) (5) 

for scattering states with real energy E = h 2 k 2 /2m = E k , which are incident from the 
left (k > 0), 

m': h ill < 6 » 



_ T 2 e i i x 

These states are not orthogonal, in spite of the reality of E, because the Hamiltonian 
H c is not Hermitian. The wavenumber q obeys 

E + ihj/2 = ti 2 q 2 /2m, (7) 

with Img > for boundedness, while R12 and T 12 are reflection and transmission 
amplitudes yet to be determined for the ground and excited state channels. 

To obtain the form of for < x < L we denote by |A + ) and |A_) the 

(unnormalized and nonorthogonal) eigenstates of the matrix ~ (q J] 7 ) corresponding to 
the eigenvalues X±, 



1 % 



' A± > - U/n) ' < 9 » 

(Formally we assume that A + 7^ A_. The case A + = A_ may be treated by taking the 
limit.) For < x < L, one can write as a superposition of |A±) in the form 

V2n& k (x) = C ++ \X + )e lk+x + C_+|A_)e ifc - x (10) 
+ C + ^\X + )e- ik+x + C—\X^)e- ik - x , < x < L, (11) 

which, at variance with the semi-infinite laser case must contain now growing 
exponentials in addition to decaying ones. From the eigenvalue equation H c <&k = E^^, 
together with Ek = h 2 k 2 /2m, there follows 

k\ = k 2 - 2mX±/h = k 2 + ijm/2h =p i^ff - 4Q 2 m/2h, (12) 

with lmk± > 0. The continuity of ^(x) and its derivative at x = and x = L leads to 
eight equations with eight unknowns and to explicit but rather lengthy expressions for 
the coefficients of the wave function in the laser region, see (|l()jl . and for the transmission 
and reflection amplitudes. We shall not display them here but discuss different limits, 
approximations and numerical examples. 
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An important simplification occurs for kinetic energies above Qh/2 or, equivalently, 
for de Broglie wavelengths 2ir/k smaller than the Rabi- wavelength 2tcv/Q: the stationary 
wave may be very well approximated according to the classical idea that the atom moves 
along an unperturbed classical trajectory, without reflection, and with its internal level 
populations oscillating quantally in the laser-illuminated region. Specifically, *&(x) in 
the laser region (0 < x < L) can be approximated as the product of the internal wave 
function for the atom at rest, assuming a pure ground state at an entrance instant t — 0, 
and making the substitution t = x/v, times the translational plane wave e tkx / '{2-k) 1 / 2 , 

e - e -W- J 008^-774] + -^L— sin[£ - 7 '/4] 

<2{X) = -—p: < V ' (13) 

W ~ in sinff JQ 2 -j 2 /4] 

whereas to the left of the laser R\ w i?2 ~ 0, and to the right the transmission amplitudes 
Ti_2 are obtained by continuity from (jUJ) and f|13j) . In particular, 

L 



T = ^k- q) L e - lL /A V sin 



(14) 



(Q2 _ 7 2/4)l/2 

The validity of this semiclassical approximation for the translational part of the 
stationary waves does not mean that every wave packet with energy components well 
above Qh behaves classically Two or more different stationary components may add up 
coherently in the time dependent wavefunction, leading to very non-classical interference 
effects, as we shall see in an example below. 

If ip{k) denotes the wavenumber amplitude that the wave packet would have as a 
freely moving packet at t — 0, then 

poo _ „ 

*(x,t)= dktfj(k)<f> k (x)e- ihk t/2m (15) 



describes the "conditional" time development of the state for an undetected atom which 
in the remote past comes in from the left in the ground state. 



3. Operation regimes 

Figure ^ shows the total detection probability, or "absorption" due to the non-hermitian 
potential, A = 1 - [T^ 2 - \R l \ 2 1 versus the incident velocity v and Rabi frequency £1 
for the stationary wavevectors A may also be understood as the probability that 
the particle be reflected or transmitted in state |2), or be detected through a photon 
emission in the laser region. We may distinguish several regions in the v—£l plane (these 
regions are also portraid in figure EJ) depending on three basic criteria, each associated 
with a dimensionless parameter: 

• Reflection: Except for very weak driving, where it is not significant, reflection is 
important for kinetic energies below or around E = Qh/2, (dashed line in figure |2J), 
but vanishes when E » Qh/2. 

• Laser intensity: Strong driving for Q/'y > 1, and weak driving for Q/'j < 1. 



Time- of- arrival measurements 



6 




• Laser beam width: Let I be the "penetration length" of the stationary wave in the 
semi-infinite laser region, that we may estimate as five times the detection delay (7| 
multiplied by the atomic velocity v, 

l = 5, g + . (16) 

For L > I the laser behaves effectively as a semi-infinite one. This regime, to the 
left of the curved solid line of figure corresponds essentially to the one examined 
in the previous paper. The exact expressions for wave functions, reflection and 
transmission amplitudes obtained in [7j are good approximations for the finite width 
laser if L > I but fail otherwise. Figure El shows the critical "temperature" versus 
L (for Cs atoms and 7 = 33.3 x 10 6 s _1 ) above which the finite width of the laser 
has to be taken into account. 

For a constant Q, at moderate to strong driving, A increases very sharply from 
zero to one at low v. This increase is not visible in figure Q because of the scale, 
but can be seen in the close-up of figure EJ The small absorption near v = is due 
to reflection into the ground state. After the rapid increase of A at low v along a 
constant Q line, there follows a plateau of perfect detection, still within the "semi- 
infinite laser" regime, see again figure H Gradually, for higher and higher velocity v and 
in strong driving conditions, the plateau is substituted by an oscillating pattern. The 
absorption surface to the right of the plateau resembles a fan, which can be explained 
in terms of the semiclassical picture and the finiteness of the laser "barrier" . The ridges 
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Figure 2. il-v plane with separations of the different regions, L/l = 1 (curved solid 
line), 2E /Ml = 1 (dashed line), and maximum absorption lines up to n — 20 (straight 
solid lines), see lfT%|) . L — 5 /mi as in figure H 



correspond to combinations of parameters where the Rabi oscillation leaves the excited 
state population at a maximum on the laser edge L, so that the decay takes place 
behind the laser, at x > L, with a detection delay I/7 from the instant when the 
classical trajectory leaves the laser. 

Each ridge corresponds to a number n+ 1/2 of Rabi oscillations between x = and 
x = L, and will be indexed by n = 0, 1, 2.... Neglecting the effect of 7 within the laser 
region, the velocity at the ridge for a given Q is given by 

m , . 

Vn = (2^TTF ( } 

Equivalently, each absorption ridge may be characterized by the straight line 

n. = (is) 

in the Q—v plane. The first one, for n — 0, corresponds to a spatial version of a 7r 
pulse, i.e., the velocity, laser-beam width, and Rabi frequency are just the right ones to 
pump the (semiclassical) atom, from the ground state at x = to the excited state at 
x = L, in half Rabi oscillation. Note from (jl7j) that the fan structure "folds upwards" 
by decreasing L, because smaller lengths require higher laser intensities to achieve full 
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Figure 3. "Critical Temperature", calculated as mvj^/K, versus laser-beam width L 
for Cesium atoms and strong driving conditions, vl = 2/y/10, i.e., the velocity for 
which the penetration length is equal to the laser beam width, I = L. 7 = 33.3 x 10 6 
s _1 , and K is Boltzmann's constant. Above this temperature the effect of the finite 
width of the laser beam cannot be ignored. 



detection. The minima at the valleys correspond to an entire number of Rabi oscillations. 
These minima are not exactly zero due to the absorption that takes place in the laser 
region. As the velocity v increases, however, this absorption decreases, as may be seen in 
figure El Finally, for velocities greater than t> , the atom is too fast for being completely 
excited during the crossing of the laser-illuminated region, and the absorption decreases 
monotonously. The solid line in figure El shows a cut of the absorption surface along a 
constant-fi line for strong driving conditions. 

The velocity window for nearly complete detection on a ridge along a constant-^ 
line is larger for smaller n. This is also clearly illustrated in figure El A simple formula 
for the window width Ayy around the maximum where absorption is better than 99% is 

1 2 



Aw < 



m 
To 



4v r . 



(19) 



yr(2n + l)J 10 107r(2n + l)' 

Notice that even in the best case, n = 0, a wave packet located essentially within one 
of these windows of full detection will have a small width Av compared to its average 
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Figure 5. Absorption versus v for £1 = 57 (solid line), and Q = 7/2. L = 5/oti. 
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velocity. In particular, for A w = 8Av, 

Av « ^— . (20) 

(2n+l)60 V ; 

An interesting feature of the transmission amplitude T 2 at the detection maxima is its 
alternating sign, —i for n even, and % for n odd, see ()14|) . We shall comment on physical 
consequence below. 

The above description is not applicable to weak driving conditions. Reflection in 
particular is less important at low velocities than for strong driving, and it vanishes for 
any v 7^ when — > 00. Also, the plateau of full detection is much narrower than for 
strong driving, and for higher velocity v along a constant-^ line, the laser does not have 
enough intensity to produce Rabi oscillations, so the absorption decays monotonously 
with v, see figures EHH and El 

The detection delay for a given incident wavenumber will in general have 
contributions from the laser region and from the region outside. The delays in the 
laser region are 2/7 for strong driving or j/fl 2 for weak driving when L > I [7], or 
fractions of these quantities for smaller laser-beam widths, whereas the detection delay 
once the (semiclassical) excited atom arrives at L is I/7. This later contribution would 
be very small for large (ideally infinite) 7, but in this limit, with all other parameters 
constant, there would not be excitation at all. One might then imagine the case where 
both 7 and Q are infinite, but with a constant ratio fi/7. In this case, however, the 
penetration length goes to zero, and the semi-infinite laser result, R\ — > —1, applies, so 
that all atoms would be reflected without being excited [7j. Another not very helpful 
limit is L —>■ 0, with all other parameters fixed, because all atoms are transmitted 
in the ground state without being detected^. "Good" direct (no delay, full detection) 
measurements are possible, but only within specific parameter ranges. For moderate 
velocities, i.e., provided L > I for the wave packet momentum components, and above 
the reflection region, the measurement on the full-detection plateau can be very efficient, 
as illustrated in Ref. [IJ. This requires, for strong or weak driving, 
Ti 1 7 2 L 

B«rf«^«5? At > < 21 > 
where At is the time span of the wave packet passage. Even when these conditions 
are fulfilled, the observed II is only approximately equal to J. Unfortunately, the 
deconvolution used in Ref. jxj to get the flux J exactly for large 7 may be impossible 
to implement with a finite-width laser at large velocities, since the penetration length 
I — ► 00 as 7 — > 00. The following section provides a way out which, again, is not 
universal but depends in fact on rather restrictive conditions. 

I One might of course insist on using the information of the very few detected atoms and normalize, 
but the result is strongly biased in favour of the slowest energy components. 
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4. An ideal distribution 

Full detection is achieved sufficiently close to the ridges of the fan structure for moderate 
to strong driving. The first maximum, corresponding to n = 0, is particularly suitable 
because it provides the largest momentum window of nearly complete absorption. Let 
us assume that Q and L may be adjusted so that most of a wave packet lies within 
an absorption maximum. As in the weak driving case of the semi-infinite laser, full 
detection comes with a price: the detection delay due to the time necessary to de-excite 
the atom, now outside the laser illuminated region. We have pointed out already that a 
short lifetime (a large 7) "solves" this problem but also creates a new one: the absence 
of excitation and therefore of detection. This suggests that the useful limit may be 
quite the opposite, namely, 7^0. The long delay can then be substracted using a 
convolution formula to obtain an "ideal" distribution that coincides with the flux when 
7^0, or, quite independently of 7, for narrow wave packets in momentum space. 

To define from the "experimental" U.(t) an idealized arrival-time distribution, we 
shall assume 

h 1 <L 1 

« <_ « 22 

hi \l v 7 

so that we can use the semiclassical and strong driving approximations neglecting the 
effect of 7 between and L, since the only significant absorption (i.e., detection) occurs 
for x > L. 

n id (t) is then defined by the convolution formula 

n = n id * w, (23) 

where 

W{t) = 7e" 7 * (24) 

is the probability density to detect a photon at time t if the atom is excited at t = in 
the laser-free region. Even though (J23|) is in fact the equation that defines iTy, there is of 
course a "classical" physical motivation for the convolution structure. If the arrival time 
of the atom, and the time of photon emission from the arrival instant were independent 
random variables, then the distribution for the sum of these two quantities would have 
precisely that form. It turns out, as discussed in Ref. [Zj that n id may become negative 
so, generically, this simple hypothesis does not hold true quantum mechanically. This 
does not invalidate Il id , but should make us cautious about its interpretation: it plays 
the role of a time of arrival density in quantum mechanics, but does not share all the 
properties of the corresponding classical density, most prominently positivity. This 
resembles the status of the Wigner function, which plays the role of a classical phase- 
space probability density but can also be negative. 

The ideal distribution Ilid is obtained from the relation between Fourier transforms, 
n id = U/W where U(u) = J dte' wt U(t) etc. From (J23) one finds 

W{v) = , (25) 

v ' w + 7 
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so that, in the time domain, 

n id (t) = n(t) + -n'(t). (26) 

7 

If the wave packet is centered at one of the absorption maxima, and the momentum 
width satisfies (fT9*j). the atoms are all pumped to the excited state, \T 2 \ ~ 1, and we 
may write 

U(t) = — l°° dx folk [ dk' ^(k)ij*(k')e' tiE ~ E ' )t/h e iiq - q ' ,)x . (27) 
2tt Jl J J 

Performing the x integral, and approximating the resulting denominator according to 

/* , , / 27m k + k' , , 

q -f„ k + (28) 



one finally gets 



n id (t) = — J J dkdk'e- l{E - E ' )t/h e i{k - k ' )L ^(k)^*(k') 
i~1 + (k 2 - k' 2 )h/2m 



(29) 



The term in square brackets tends to the kernel of the integral that gives the current 
density J, 

<*±*2» (30) 
2m ' v ; 

in several limits. Using k! = k + 5, this occurs in particular for 7 — > 0, and also for 

S/k — > independently of the value of 7, as long as is satisfied. In this later 

case, namely for a small momentum width in comparison to the momentum itself, (J3U|) 

coincides also with the kernel of Kijowski's axiomatic TOA distribution, 

(W (31) 

m 

since ~ k' and the arithmetic and geometric means in (JHUj) and (|ST|) coincide. The 
condition of small momentum-width also implies that the effect of the derivative in 
is negligible, and that, even if the last inequality of fails, J can be obtained very 
accurately by normalizing the observed results. An example is shown in figure where 
two Gaussians on the top of a single ridge have been added coherently, 

V> = + ^b, (32) 

to produce an interference pattern in the flux. 

According to (J2(J|) . the condition Av/(v) << 1 is necessarily fulfilled when the 
wave packet momenta are within one ridge of maximum detection, so that the flux and 
Kijowski's distribution will coincide, i.e., no negative values of J may be found in this 
case. In order to measure backflow the wave packet must be located within the full 
detection plateau and then proceed as in Ref. [Zj but the finite width of the laser beam 
imposes a maximum speed for this method to work. Another route would be to use 
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Figure 6. Flux (solid line), Ilid (dashed line), and normalized (dots), for 
the coherent sum of two Gaussians with average velocities v± = 167.05m/s and 
Vi — fi + .9/Ltm/s, and Ax = 4233/im. They become minimal-uncertainty-product 
states when their centers arrive at L = 5/xm. 51 = 104.43 x 10 6 s _1 . 



different ridges of maximum detection for each gaussian. The two gaussians should be 
on odd or even n ridges but not on contiguous ones because of the sign change of the 
transmission amplitude mentioned below (|20j) . 



5. Conclusions 



In this paper we have analyzed the measurement of the arrival time of an atom at a given 
location by means of a laser that illuminates a finite region. Repeating the experiment 
many times for a given atomic preparation produces an "observed distribution" of first 
fluorescence photons II, that will in general be distorted with respect to the "axiomatic" 
distribution of Kijowski 11^ or the flux J (both evaluated for the isolated atom, i.e., 
without the laser) because of detection delay, and atomic reflection or transmission 
without photon detection. (In general even if the three problems are made negligible 
II is only approximately equal to these "ideal" quantities.) It is possible to correct for 
the distortions due to detection delays, at least in certain limits, and obtain the current 
density by deconvolution. In the present model the laser beam intensity has sharp edges 
to facilitate the analytical examination but very similar results have been obtained for 
a more realistic Gaussian profile, see figure [7| 
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Figure 7. Absorption versus velocity v and fio for L = 5fim and position dependent 

Rabi frequency, fl(x) — LVl^e 5P /5(2n) ' ; S = 0.529/im, xq = 2.5/im, and 
7 = 33.3 x 10 6 s _1 . The calculation is performed with a transfer matrix method. 
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